Hindered rolling and friction anisotropy in supported carbon nanotubes 
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Carbon nanotubes (CNTs) are well known for their exceptional thermal, mechanical and electrical 
properties. For many CNT applications it is of the foremost importance to know their frictional 
properties. However, very little is known about the frictional forces between an individual nanotube 
and a substrate or tip. Here, we present a combined theoretical and experimental study of the 
frictional forces encountered by a nanosize tip sliding on top of a supported multiwall CNT along 
a direction parallel or transverse to the CNT axis. Surprisingly, we find a higher friction coefficient 
in the transverse direction compared with the parallel direction. This behaviour is explained by a 
simulation showing that transverse friction elicits a soft hindered rolling of the tube and a frictional 
dissipation that is absent, or partially absent for chiral CNTs, when the tip slides parallel to the 
CNT axis. Our findings can help in developing better strategies for large-scale CNT assembling and 
sorting on a surface. 



Carbon nanotubes (CNTs) are well known for their ex- 
ceptional thermal, mechanical and electrical properties 
[THE]. For many CNT applications it is of the foremost 
importance to know their frictional properties. However, 
very little is known about the frictional forces between 
an individual nanotube and a substrate or tip. Here, we 
present a combined theoretical and experimental study of 
the frictional forces encountered by a nanosize tip slid- 
ing on top of a supported multiwall CNT along a direc- 
tion parallel or transverse to the CNT axis. Surprisingly, 
we find a higher friction coefficient in the transverse di- 
rection compared with the parallel direction. This be- 
haviour is explained by a simulation showing that trans- 
verse friction elicits a soft hindered rolling of the tube 
and a frictional dissipation that is absent, or partially 
absent for chiral CNTs, when the tip slides parallel to 
the CNT axis. Our findings can help in developing bet- 
ter strategies for large-scale CNT assembling and sorting 
on a surface. 

Frictional forces are responsible for an enormous va- 
riety of phenomena in our everyday life including play- 
ing the violin, dancing the tango and driving a car [7J. 
The same wide spectrum of friction applications can be 
found in nanobiotcchnology. In particular, the frictional 
forces between concentric CNTs have attracted consid- 
erable interest for CNT potential use in nanomechanical 
[H [3] and nanoelectrical systems [HI H] • Experiments and 
theory agree that the translational and rotational inter- 
layer forces between suspended concentric nanotubes are 
extremely small, with values of the translational shear 
strength ranging between 0.2 and 0.04 MPa (refs 10, 11) 
and values of the rotational shear strength varying from 
0.85 to 0.05 MPa (refs 12, 13). The frictional behaviour 
of a nanotube deposited on a surface can clearly be very 



different. Further effects need to be considered, such as 
the possibility of rolling and the role of the adhesion force 
between the nanotube and the substrate. Atomic force 
microscopy (AFM) is an excellent tool to manipulate ob- 
jects and investigate forces at the nanoscale. Previous 
AFM studies in air have investigated the sliding of a 
CNT on a graphite surface by CNT end-on pushing with 
an AFM tip [U]. More recently, a CNT tip attached 
to a conventional AFM probe was slid along a nanotube 
suspended over a trench [IS] . Both studies found for the 
carbon-carbon interface a shear stress of a few mcgapas- 
cals. No friction studies are available about the role of 
CNT radius, CNT structural properties, sliding velocity 
and especially sliding direction relative to the tube axis. 
As CNTs are made of concentrically rolled-up graphene 
sheets, they might be expected to yield an essentially 
isotropic nanofriction when they slide on a surface par- 
allel (longitudinal sliding) or perpendicularly (transverse 
sliding) to their axis. 

Here, we present a joint experimental and theoreti- 
cal investigation of the frictional forces encountered by 
a nanosize tip sliding on top of a supported multiwall 
CNT along a direction parallel or transverse to its axis, 
as a function of normal load, sliding velocity and tube 
radius. Experimentally, the longitudinal shear strength 
is found to remain almost constant with varying CNT 
radius, whereas the transverse shear strength decreases 
with increasing radius, reaching the value of the longitu- 
dinal shear strength for radii larger than 10 nm. Experi- 
ments and theory find that for a fixed total normal force, 
that is, applied normal load plus adhesion force, the fric- 
tion force, Fp, in the transverse direction is higher than 
in the longitudinal direction. In particular, the theo- 
retical simulations find that the transverse-longitudinal 
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FIG. 1: Friction measurements on CNTs. a, Topography image of a CNT of 7 nm radius. The fast scanning direction of 
the AFM tip is indicated by an arrow (X direction). b,c, Friction images of the highlighted longitudinal (b) and transverse (c) 
sections of the nanotube. d,e, Topography and friction force profile across the CNT. The topography profile (black solid line) 
is along the white solid lines in b,c. The friction force profile (red solid line) is the average profile inside the area delimited by 
the dotted line in b,c. The friction force profile in d is along the Y direction and the one in e is along the X direction, fh, 
Topography images of three other CNTs of radius 11, 9 and 6 nm, respectively, where friction measurements were carried out 
inside the areas delimited by solid lines. 



friction anisotropy, a = Fp an /F F 9 , can be as large 
as 20 in non-chiral CNTs and about 2-3 in chiral CNTs. 
The experiments, carried out on multiwall CNTs with 
random chirality, show that a can be as large as 3. This 
unanticipated behaviour is explained by the theoretical 
simulations. We find that transverse sliding probes not 
only the elastic deformations of the stiff carbon-carbon 
bonds of the graphitic basal planes, common to the lon- 



gitudinal sliding, but also a much softer overall swaying 
motion of the tube, akin to a hindered rolling. This added 
softness is the source of the extra transverse tip dissipa- 
tion and increased friction. 

The multiwall CNTs used in this study are produced 
by a catalytic chemical vapour deposition method and de- 
posited on a flat silicon substrate. The structural proper- 
ties of the CNTs were characterized by transmission elec- 
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FIG. 2: Frictional forces for transverse and longitudinal 
sliding. a,b, Frictional force as a function of the normal load 
for silicon and for transverse (a) and longitudinal (b) sliding 
on top of a CNT, measured as described in Fig. 1. The CNT 
radius is ~ 3.5 nm in the longitudinal section and ~ 5 nm in 
the transverse section. The sliding velocity was 0.8 fim/s. The 
solid lines are fits to the data by using Ff = {i{Fn + F a dhY ■ 
The error of Ff is obtained from the fluctuations in the values 
of friction in the investigated areas on top of the nanotubes 
(see Supplementary Information for more details). 



tron microscopy, which indicated that these CNTs have 
several defects and a constant ratio of external to inter- 
nal radii of RNr/Rmt = 2.8 ± 0.7 (see Supplementary 
Fig. S4). The topography and the friction properties of 
the CNTs were characterized by AFM. To measure the 
transverse and longitudinal friction force on a nanotube, 
we looked for a statistical number of CNTs on the silicon 
surface and we took simultaneous topographical and fric- 
tional images in different locations of the nanotube and 
at different scales. In Fig. 1, we show some examples 
of CNTs that have been measured. We usually look for 
nanotubes that lie on the substrate with some sections 
parallel and some other sections perpendicular to the tip 
scan direction. We then zoom into the desired section 
and we acquire a friction and topography image (Fig. 



lb,c). This method also has the advantage that for each 
friction/topography measurement, we acquire during the 
same scan FF of both the nanotube and the silicon sub- 
strate, which can be used as our reference system. Typ- 
ical cross-sections of the friction and topography images 
for longitudinal and transverse measurements are shown 
in Fig. ld,e. Figure 2a,b shows a typical plot of Fp 
acquired on top of a CNT as a function of the normal 
load, Fn, for transverse (Fig. 2a) and longitudinal (Fig. 
2b) sliding. In each graph, we also add the results for 
the silicon substrate, which is the reference surface. As 
a first approximation, Fp between the AFM tip and the 
silicon substrate or between the tip and the CNT, could 
be expressed as: 



Fp — /U(Fjv + F a dh) 



(1) 



where \x is the friction coefficient and F a dh is the adhe- 
sion force between the AFM tip and the investigated sam- 
ple. However, it is well known that this linear relation- 
ship becomes a quasi-linear relationship at the nanoscale, 
where the "single contact" geometry is very likely. In 
fact, the friction force is proportional to the tip-sample 
contact area, A(Fjq + F a dh)' 



F F = a- A(F N + F adh ) 



(2) 



where the proportionality factor is called shear 
strength, a (ref. 16). The contact area between a spher- 
ical tip (AFM tip) and a cylinder (CNT) or a flat sur- 
face (silicon substrate) can be expressed as a function 
of Fjsr, F a dh, elastic moduli, tip radius, Ru P , and cylin- 
der radius, Rnt, by using continuum mechanics theo- 
ries, such as the Hertz theory [17] . The details of the 
calculation of A(F^ + F a dh) are reported in the Sup- 
plementary Information. As a result, one obtains that 
Ff = £i(Fn + Fadh) 2 / 3, ■ Depending on the nature of the 
adhesive forces, this equation can be slightly different; 
however, a discussion of the different contact mechan- 
ics models is beyond the scope of this letter and we re- 
fer the readers to refs 18 and 19 for details. Figure 2 
seems indeed to indicate a 2/3 power-law dependence, 
but more importantly it clearly shows that the frictional 
forces are very different when the tip slides parallel or 
perpendicularly to the CNT long axis. The slope of the 
Fp versus Fn curve for the transverse CNT is similar to 
the slope measured on Si and it is about three times the 
slope measured on the longitudinal CNT. On the other 
hand, F ad h remains almost the same for the two sliding 
directions (5.2 nN for longitudinal sliding and 3.3 nN for 
transverse sliding, as obtained by fitting the data with 
Ff = fi(FN + F a dh) 2 ^ 3 '■ No substantial differences have 
been detected in the Fp versus Fn curves acquired on 
the same CNT and for the same sliding direction but for 
different sliding velocities, in the range v — 0.8 — 4 /im/s. 

To further analyse the experimental data and compare 
the friction behaviour of CNTs with different radii, we 



4 



a 2.5 



2.0 



S 15 



1.0 



0.5 



i: 



# Transverse 
I Longitudinal 



6 8 10 12 14 

Nanotube radius (nm) 



b 15 



# Transverse 
■ Longitudinal 



i 



i 



8 10 
Nanotube radius (nm) 



12 



1-1 



FIG. 3: Shear strength and adhesion force. a,b, Shear 
strength (a) and adhesion force (b) for transverse and longi- 
tudinal sliding on top of a CNT as a function of the tube ex- 
ternal radius. The sliding velocity was 2 fj,m/s for all of these 
measurements. The error on the shear strength is determined 
by fitting FF versus FN with equation (2), as described in the 
Supplementary Information. The error on the adhesion force 
is obtained by the fit of Ff versus Fn with equation (1). 



derive from the experimental Fp versus Fn curves the 
shear strength, as introduced in equation (2) and F adh 
(see Supplementary Information and the Methods sec- 
tion). Figure 3a,b shows a and F ad h, respectively, as a 
function of the CNT radius for transverse and longitudi- 
nal sliding, at fixed velocity (2 /im/s). The longitudinal 
shear strength (~0.5 GPa) is independent of the tubes' 
radii in the range between 4 and 12 nm. This value is 
three/four orders of magnitude larger than a measured 
between the walls of the CNTs during telescope-pulling 
experiments. However, it is only two orders of magni- 
tude larger than the values (about 4 MPa) found in ex- 
periments where a CNT slides on graphite or a CNT 
tip slides on another CNT in air. Furthermore, from 
our experiments carried out on highly ordered pyrolytic 



graphite in air with a Si tip we find ohopg = 3 — 60 
MPa, in good agreement with previous results [20] . The 
presence of water, the tip-CNT contact geometry [T5] 
and the presence of defects in the CNT seem thus to be 
the reason for a large value of er measured during our 
experiments on CNTs. By comparing longitudinal and 
transverse sliding, the experiments indicate that in the 
transverse direction a is about three times larger than in 
the longitudinal direction for a tube radius of 4 nm. For 
larger tubes, a slightly decreases in the transverse direc- 
tion, with some fluctuations that could be related to a 
different CNT chirality and/or to the interaction of the 
CNT with the substrate, as suggested by the values of a 
for the radius of 11 nm for which a = 0.825 GPa when 
F adh = 11.2 nN and a = 1.781 GPa when F adh = 5.6 nN. 

To understand the observed frictional behaviour and 
the surprising transverse-longitudinal anisotropy, we con- 
ducted extensive molecular-dynamics simulations. Sev- 
eral model multiwall CNTs were built, with both arm- 
chair and chiral geometries and deposited on ideally flat 
crystalline substrates, the role of which is simply to sta- 
bilize the CNTs. The results presented here are for two 
different four-wall nanotubes, one fully armchair (26, 26)- 
(31,31)-(32, 40)-(41, 41) non-chiral tube and one (4, 46)- 
(32, 32)-(37, 37)-(18, 63) chiral tube. Well-tested force 
models (see the Methods section) were found to stabilize 
tubes with reasonable geometries and compliances (Fig. 
4a, b). Each tube was approached by a nanosized dia- 
mond tip, controlled by a cantilever with three spring 
constants, exerting a variable load Fjq in the nanonew- 
ton range and laterally sliding back and forth with a 
speed v ~ 2 m/s along longitudinal and transverse di- 
rections relative to the tube axis. As in the experi- 
ments, the sliding amplitude was of the order of one 
nanometre and the frictional dissipation was extracted 
from the lateral force/displacement parallelogram area. 
The simulations clearly showed that the tubes deformed 
in two distinct manners during the tip pressing and slid- 
ing. First, we observe the formation of a dimple when 
the tip is on top of the tube centre (Fig. 4c, h). Sec- 
ond, we find a strong shape asymmetrization accompa- 
nied by an attempted lateral rolling of the tube (blocked 
in our case to prevent artificial runaway rolling caused 
by periodic boundary conditions) when the tip is off the 
tube centre (Fig. 4a,b,e-g). During motion, the lon- 
gitudinal sliding simply displaced the dimple under the 
tip and along the tube. The transverse sliding instead 
gave rise to a further large deformation and lateral sway- 
ing of the tube. Under lateral pushing by the tip, the 
tube deformed as a whole, attempting a genuine hin- 
dered rolling back and forth in time. The calculated 
tip-nanotube frictional forces, determined by the dissi- 
pated energy over many cycles, are shown in Fig. 5, for 
both non-chiral and chiral nanotubes and for longitudinal 
and transverse motion. As observed during the experi- 
ments, the simulations show large transverse-longitudinal 
friction anisotropy. The larger transverse dissipation is 
clearly shown by simulation to be connected with the 
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FIG. 4: Molecular dynamics simulation of the tipnanotube interaction. a,b, Non-chiral tube without tip (a) and under 
2.5 nN normal force (b). c, The sectional view of the dimple, d, View of the outer tube (armchair) . eg, Chiral tube without tip 
(e) and under 1.9 nN (f) and 2.5 nN (g) with the tip pressing on two different spots, h, Sectional view, under 2.5 nN. i, Outer 
view of the chiral outer tube. Smaller sizes of both CNT and tip account for numerically smaller frictional forces compared with 
the data of Fig. 2. 



hindered rolling soft motion. Hindered by substrate ad- 
hesion, the tendency to roll decreases for increasing radii, 
consistent with the experimentally observed anisotropy 
decrease at larger radii. It should be noted that there is 
a six-order-of-magnitude difference between the theoret- 
ical and the experimental sliding speeds (see Methods). 
However, in both cases the energy dissipation occurs by 
transfer of energy into the nanotube modes: smoothly in 
the simulation and by "kicks" during the probable stick 
and slip motion in the experiments. 

Interestingly, the simulations show an anisotropy 
pTran J pLong ag Y &r g e as 20 in the fully armchair non- 
chiral tubes, but it drops to about 2 when the outer 
tube is chiral. Owing to chirality, the tip causes in a 
screw-like fashion some hindered rolling even for longitu- 
dinal sliding, owing to the absence of left-right symme- 
try. This suggests that our experimental multiwall CNTs 
were generally randomly chiral |21) . It should be men- 
tioned that structural defects present in the CNT might 
also contribute in coupling the transverse and longitudi- 
nal dissipation modes in non-chiral CNTs, giving rise to 
an anisotropy coefficient smaller than 20. Dissipation due 
to motions outside the CNT and the occurrence of stick- 
slip probably also mitigates the frictional anisotropy. 

We have shown an anisotropy in the shear strength of 
CNTs in the transverse and longitudinal direction, which 
has its origin in the soft lateral distortion of the tube 
when the tip-tube contact is moving in the transverse di- 



rection. This friction anisotropy should be observable in 
very different geometries/situations, for example, when 
a CNT slides on a substrate parallel or perpendicularly 
to its axis. Our findings could help in developing better 
strategies for chirality sorting, large-scale self-assembling 
of surface-deposited nanotubes [35] , design of CNT adhe- 
sives [5] and nanotube-polymer composite materials [BJ. 
The conclusions of this letter more generally affect the 
ability to manipulate nano-objects on a surface. 



Methods 

AFM. The topography and frictional force images 
were collected simultaneously using a Veeco CP-II 
AFM with a silicon tip of about 60 nm (PPP-CONT, 
NanoWorld) in air at room temperature, at a scanning 
speed ranging between 0.8 and 4 /jm/s, for different 
normal force values. The normal force was between 8 
and 2 nN initially and then decreased until the tip was 
out of contact with the sample. The cantilever nor- 
mal spring constant of O.llN/m was calibrated using the 
method of Sader and colleagues 23.. The lateral force 
was calibrated using the wedge method on the Mikro- 
masch TGG01 silicon grating [24] (Supplementary Figs 
SI and S2). The diameter of each section of CNT was 
inferred from the height of the imaged CNT lying on the 
silicon substrate. The contact area between the AFM tip 
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FIG. 5: Molecular dynamics simulation of tipnanotube friction force versus normal load, a, Non-chiral nanotube 
(very large transverselongitudinal anisotropy). b, Chiral nanotube (more realistic transverselongitudinal anisotropy of about 2). 
The frictionload behaviour is to a good approximation linear (solid lines) and not a (2/3) power law as in the experiments, 
as the small size of the simulated tip and of the tipCNT contact makes continuum contact mechanics invalid in this limit. 
Simulated friction extrapolates to zero at zero load, because the experimental zero-load situation corresponds to a finite load in 
the theory, the difference between the two corresponding to the experimental tip pull- off force, which measures adhesion, omitted 
in the simulation. The error bars indicate the friction fluctuations between different sliding loops. 



and the silicon substrate or between the tip and a CNT 
was evaluated using the modified Hertz theory describ- 
ing, respectively, the contact between a spherical tip and 
a planar substrate, or the contact between a spherical tip 
and a cylindrical tube [THl [55] ■ 

To evaluate Fp on the very top of the nanotube, after 
imaging the friction and topography of the desired sec- 
tion of the CNT and the surrounding silicon substrate, 
we zoom in on an area at the top of the tube where the 
topography and the friction look flat. For example, for 
a CNT of radius 7 nm, we zoomed in on an area of 20 
nm x 100 nm in the centre of the tube. By considering 
the convolution of the AFM tip (radius 60 nm) with a 
CNT of 7 nm, a width of w con voiuted — ±10 nm at the 
top of the imaged tube corresponds roughly to a width 
w re ai = ±1 nm on top of the real nanotube. See Sup- 
plementary Information for more details on the evalua- 
tion of the contact area, friction measurements and error 
analysis (Supplementary Figs S5 and S6). 

Simulations. In the molecular-dynamics simulations, 



CNTs were deposited on crystalline substrates, assumed 
for definiteness to be gold (111) (in place of the real 
and more complex oxidized Si surfaces). The outer/inner 
tube radius ratios were substantially smaller in our sim- 
ulations than in the experiments, implying first of all nu- 
merically smaller frictional forces and also a greater soft- 
ness and a quantitatively larger tip-induced deformation 
from the cylindrical shape, as shown in Fig. 4. Periodic 
boundary conditions were assumed parallel and perpen- 
dicular to the tube. For the substrate-nanotube interac- 
tion, a Lennard- Jones Au-C potential j26| was used with 
a strength decreased by a factor of 10, an interaction that 
does not excessively alter the CNT cylindrical shape; all 
substrate atomic positions were assumed to be rigid, jus- 
tified by the much larger substrate rigidity relative to 
the CNTs. Standard empirical Brenner potentials [27] 
and Kolmogorov-Crespi potentials [23] were used among 
CNT atoms. On the basis of the justification that the re- 
pulsive forces are short ranged and primarily responsible 
for the normal load, the tip-tube potential was assumed 
to be equal to the repulsive part of the Lennard- Jones 



potential, V(r) = Ae{a/r) (rcf. 12), with e = 0.004 eV 
and a = 0.328 nm. A Berendsen thermostat was applied 
to the unconstrained atoms of the outer tube. The tem- 
perature was set at T = 50 K, different from the room 
temperature at which experiments were carried out, to 
avoid large fluctuations of the frictional force related to 
the small CNT sizes. Artificial massive lateral rolling of 
the tube, caused by the infinite tip replicas implied by 
the periodic boundary conditions and prevented in ex- 
periments by resistance to tip pushing by the strong sub- 
strate adhesion, was blocked in the simulation by freez- 
ing the relative carbon-substrate positions for a minimal 
number of outer tube strands in closest contact with the 
substrate (yellow coloured in Fig. 4). The rolling can be 
both leftwards and rightwards (Fig. 4f,g). The Young 
modulus ratio of silicon to CNT is roughly 170:30, mean- 
ing that the (silicon) tip and substrate deformations, al- 
though surely finite, are negligible compared with those 
of the CNT; thus, we chose a 525-atom rigid diamond 
tip to slide on the tube surface. The tip is connected to 
three supports (X, Y and Z) by springs with constants 
k x = 50 N/m along the sliding direction, k y — 500 N/m 
along the other direction in the plane of the substrate and 
k z = 2 N/m along the normal direction. The frictional 
force was calculated from the lateral force/displacement 
parallelogram area (see Supplementary Fig. S7). The 
speed of the support X was set to be \u\ = 2 m/s, sweep- 
ing forth and back, and Y and Z move together with the 
tip. This tip speed is many orders of magnitude larger 
than in the experiments, and in fact the simulated fric- 
tion is smooth without stick and slip, the latter probably 
occurring in the experiments. Despite that, we do not ex- 
pect the difference to affect the validity of our mechanism 
for the frictional anisotropy, because the CNT hardness 
against longitudinal shear and its softness against trans- 
verse shear found in the simulation should affect both 
frictional regimes in a similar manner. 
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